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Escape rate from a metastable state weakly interacting with a heat bath driven by an 
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Based on a system-reservoir model, where the reservoir is driven by an external stationary, Gaus- 
sian noise with arbitrary decaying correlation function, we study the escape rate from a metastable 
state in the energy diffusion regime. For the open system we derive the Fokker-Planck equation 
in the energy space and subsequently calculate the generalized non-Markovian escape rate from 
a metastable well in the energy diffusion domain. By considering the dynamics in a model cubic 
potential we show that the results obtained from numerical simulation are in good agreement with 
the theoretical prediction. It has been also shown numerically that the well known turnover feature 
can be restored from our model. 

PACS numbers: 05.40.-a, 02.50.Ey, 82.20.Uv 



I. INTRODUCTION 



Ever since Kramers proposed his seminal work [ij for 
chemical reaction in terms of the theory of Brownian 
motion in phase space, the model and its several vari- 
ants remain ubiquitous in many areas of natural sciences. 
Through the years it has been a subject of several the- 
oretical 1 HE II II II and experimental 0, Elll 
investigations for understanding the nature of activated 
rate processes. In the majority of these treatments, one is 
essentially concerned with a thermally equilibrated bath, 
which simulates the reaction coordinate to cross the acti- 
vation energy barrier. The inherent noise of the medium 
is of internal origin, which implies that the dissipative 
force which the system experiences in course of its mo- 
tion in the medium and the stochastic force acting on 
the system as a result of the random impact from the 
constituents of the medium arises from a common mech- 
anism. From a microscopic point of view, the system- 
reservoir Hamiltonian description [13, [l3 suggests that 
the coupling of the system and the reservoir coordinates 
determines both the noise and the dissipative terms in 
the Langevin equation describing the motion of the sys- 
tem and therefore these two entities get related through 
a fluctuation-dissipation relation |14| . which is the char- 
acteristics of a thermodynamically closed system in con- 
trast to the systems driven by external noise [IBI lla |. 
However, when the reservoir is modulated by an external 
noise, it is likely that it induces fluctuations in the polar- 
ization of the reservoir. These fluctuations in turn may 
drive the system in addition to the usual internal noise 
of the reservoir. Since the fluctuations of the reservoir 
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crucially depends on the response function; one can en- 
visage a connection between the dissipation of the system 
and the response function of the reservoir due to external 
noise, from a microscopic standpoint \u\ . At this point 
it is important to mention that, a direct driving of the 
system coordinate breaks the fluctuation-dissipation re- 
lation and can generate biased directed motion that can 
be seen in ratchets and molecular motors [l^- di the 
other hand bath modulation by an external noise agency 
preserves the fluctuation dissipation relation, as a result 
of which well known Kramers' turnover feature can be 
restored. 

In many cases involving chemical systems the Marko- 
vian representation of the Langevin equation is not valid. 
In the Markovian description the time scale associated 
with the motion of the thermal bath is much shorter than 
any relevant molecular time scale. This assumption is 
practically never realized in cases where the system co- 
ordinate is a molecular vibrational coordinate, because 
the correlation time associated with the thermal bath is 
usually much longer than a typical molecular vibrational 
period. This observation is of no consequence for the 
escape rate in the strong and moderate friction cases, 
where the particle is considered to be essentially in ther- 
mal equilibrium within the well and where the dynamics 
takes place only near the barrier top. Non-Markovian ef- 
fects may be important also for barrier crossing dynam- 
ics, however, this depends on the relation between the 
barrier frequency (renormalized by the presence of the 
friction) and the friction coefficient. But in the low fric- 
tion limit, where the well dynamics is important, energy 
accumulation becomes the rate determining step. Fur- 
thermore, reactions occurring under the nonequilibrium 
situation, the well dynamics becomes crucial and become 
dominant in the low friction regime. Obviously the well 
dynamics is governed by energy accumulation and relax- 
ation processes. In addition to Kramers' treatment in 
the low friction, there are several treatments that deal 
such situation, among which Zwanzig [13, using the as- 
sumption that reservoir is always in thermal equilibrium, 



developed a procedure for reducing the classical Hamil- 
ton's equations of motion for a one dimensional particle 
interacting with a non-Markovian heat bath. The escape 
of a particle from a potential well has been treated us- 
ing a generalized Langevin equation in the low friction 
limit by Carmeh and Nitzan (CN)|23]. Thereafter the 
detailed classical analysis reveals that the rate is signif- 
icantly modified by memory effects when compared to 
corresponding Kramers' theory. 

While nonequilibrium, nonthermal systems have also 
been investigated phcnomenologically by a number of 
workers in several contexts El HH El IM El Ei , these 
treatments concern mainly with direct driving of the sys- 
tem by an external noise or a time dependent field, e.g., 
for examining the role of color noise in stationary proba- 
bilities |2j| , pro perties of nonlinear systems ^^ , nature 
of crossover |23|, effect of monochromatic noise |2J|, to 
study the chemical reaction dynamics in an anisotropic 
solvents [25] . In the present paper we consider a system- 
reservoir model where the reservoir is modulated by an 
external noise. Our object here is to explore the role of 
reservoir response on the system dynamics and to calcu- 
late the generalized escape rate from a metastable state 
for a nonequilibrium open system in the energy diffusion 
regime. 

A number of different situations depicting the modu- 
lation of the bath may be of physically relevant. Though 
the dynamics of a Brownian particle in a uniform solvent 
is well-known, it is less clear when the response of the 
solvent be time dependent, as in the case of the dynam- 
ical properties of a suspension in a liquid crystal when 
projected on to an anisotropic stochastic equations of 
motion, or in the diffusion and reaction in supercriti- 
cal liquids and growth in living polymerization [23, Eg] . 
Also the space dependent friction may be realized from 
the presence of a stochastic potential in the Langevin 
equation [23. As another example, we consider a sim- 
ple unimolecular conversion from A^B, say an isomer- 
ization reaction. The reaction can be carried out in a 
photochcmically active solvent under the influence the 
external fluctuating light intensity. Since the fluctuation 
in the light intensity results in the fluctuations in the po- 
larization of the solvent molecules, the effective reaction 
field around the reactant system gets modified [Ig. In 
passing we mention that the escape rate in the energy 
diffusion regime is just not a theoretical issue today but 
has been a subject of experimental investigation over the 
last two decades |ll|. 

The outlay of the paper is as follows. In section II 
we discuss a system-reservoir model where the latter is 
modulated by an external noise and establish an impor- 
tant connection between the dissipation of the system 
and the response function of the reservoir due to the ex- 
ternal noise. The stochastic motion in energy space and 
Fokker-Planck equation has been constructed in section 
III. We solve the problem of energy diffusion controlled 
rate processes in section IV. An explicit example with a 
cubic potential is worked out to illustrate the theory in 



section V. The paper is concluded in section VI. 



II. THE MODEL: HEAT BATH MODULATED 
BY EXTERNAL NOISE 

We consider a classical particle of mass M bilinearly 
coupled to a heat bath consisting of N harmonic oscilla- 
tors driven by an external noise. The total Hamiltonian 
of such a composite system can be written as [13, Il3| 



2 1 ^ f 2 ^ 

Z— 1 ^ * '^ 

(2.1) 
In the above equation, x and p are the co-ordinate and 
momentum of the system particle, respectively and V(x) 
is the potential energy of the system. {qi,Pi) are the 
variables for the i-th oscillator having frequency w^ and 
mass rrii. gi is the coupling constant for system-bath 
interaction. Hint is the interaction term between the 
heat bath and the external noise, t{t) with the following 
form 



1 ^ 
Hint = T^y^Kjgieft). 



(2.2) 



The type of interaction we have considered between the 
heat bath and the external noise. Hint is commonly 
known as the dipole interaction |29l |. In Ea. H2.2(l Ki de- 
notes the strength of interaction. We consider e{t) to be 
a stationary, Gaussian noise processes with zero mean 
and arbitrary correlation function 



(e(t))e = and {e{t)e{t'))e = 2D^{t - t') 



(2.3) 



where D is the external noise strength, '^{t — t') is the 
external noise kernel and (. . .)e implies the averaging over 
the external noise processes. 

Eliminating the bath degrees of freedom in the usual 
way (and putting M and mi equal to one) we obtain the 
following generalized Langevin equation 



X = v^ 

V = -^- f dt'^{t-t')v{t') + f{t)+i:{t){2A) 



dx 







where 



N 

l{i)^^9'i^'icoBujit 
1=1 



(2.5) 



and f{t) is the thermal fluctuation generated due to 
system-reservoir interaction and is given by 



N 



f{t) = ^.9i{te(0) - giX {Q)]ujf COS ujit + Vi{0)uji sin u!it}. 
1=1 

(2.6) 
in Ea. H2.4|l . 7r(t) is the fluctuating force term generated 
due to the external stochastic driving e(i) and is given 
by 



7r(i) = - / <^(t - t')e{t')dt', 
Jo 



(2.7) 



where 



^{t) 



N 

y^^giKjUJiSmujit. 



(2.8) 



The form of Eq. H2.4|l therefore suggests that the system 
is driven by two forcing functions /(t) and 7r(i). The 
initial conditions of the bath oscillators for a fixed choice 
of the initial condition of the system degrees of freedom 
determines f{t). To define the statistical properties of 
f{t), we assume that the initial distribution is one in 
which the bath is equilibrated at i = in the presence of 
the system but in the absence of the external noise agency 
such that (/(t)) = and (/(i)/(i')> = keT-fit - t'). 

Now, at t = 0-1-, the external noise agency is switched 
on and the bath is modulated by e{t). The system is 
governed by Eq. (|2.4I) , where apart from the internal noise 
f{t), another fluctuating force 7r(i) appears, that depends 
on the external noise e{t). Therefore, one can define an 
effective noise f(0[= f{t) + vr(i)] whose correlation is 
given by 



density of the bath frequencies, the coupling of the sys- 
tem with the bath, the coupling of the bath with the ex- 
ternal noise, and the external noise itself. Equation 1)2. 7|) 
is the reminiscent of the familiar linear relation between 
the polarization and the external field, where tt and e 
play the role of the former and the latter, respectively. 
ip{t) can then be interpreted as a response function of the 
reservoir due to external noise e{t). The very structure 
of n{t) suggests that this forcing function, although orig- 
inating from an external force, is different from a direct 
driving force acting on the system. The distinction lies at 
the very nature of the bath characteristics (rather than 
system characteristics) as reflected in the relations (|2.7|) 
and ifTH)) . 

In order to obtain a finite result in the continuum limit, 
the coupling functions gi — g{Lo) and Ki = k{uj) are cho- 
sen |30| as g{u!) — go/ y^ut and k{uj) — KQUty^. Conse- 
quently 7(i) and ip{t) reduce to the following forms: 



mmt'))) - kBTj{t-t') + 2D f dt" f dt"'^{t^t") 

Jo Jo 

(2.9) 



x^(t'-t"')*(t"-t"') 
C{t ~ t') (say) , 



along with ((^(i))) — 0, where ((. . .)) means we have 
taken two averages independently. While writing (|2.9|) 
we made the assumption {{^{t)^{t'))) ~ C{t — t') which 
cannot be proved unless the structure of Lp{t) is explicitly 
given. However, as we shall see in section V and in the 
following discussion that it is a valid assumption (see, 
Eqs.(2.15-2.17) and (5.5)) for a particular choice of the 
coupling coefficients g{u)) and n{Ld) (see Eqs.(2.10) and 
(2.11)) and for a stationary external noise processes (see 
Eq.(5.1)). It should be emphasized that the above rela- 
tion 1)2. 9|l is not a fluctuation-dissipation relation due to 
the appearance of the external noise intensity. Rather it 
serves as a thermodynamic consistency condition. 

Let us now digress a little bit about 7r(i). The statisti- 
cal properties of 7r(i) are determined by the normal-mode 



7(i) = ^ I duj^{uj) cos ujt 



(2.10) 



and 



ip{t) ~ goK-o / did &{uj)uj sin ujt, (2-11) 

where go and Kq are constants and 1/tc is the cutoff fre- 
quency of the oscillator (jc may be characterized as the 
correlation time of the bath [Ij| and for Tj, ^- we obtain 
(5-correlated noise process). &{uj) is the density of modes 
of the heat bath which is assumed to be a Lorentzian 



&{lu) 



TTTc{uj^ 



(2.12) 



This assumption resembles broadly the behavior of the 
hydrodynamical modes in a macroscopic system |3l| . 
This form of density of modes, along with the expressions 
of g{uj) and k{uj), allows us to write for the expression of 
ip{t) as 



(^(t) = (goKo/rc)exp(-t/Tc). (2-13) 

From Ea. (|TTn|l and Ea. l|TTT|) one obtains [13 



di go 1 .,^ 

— = ip{t). 

dt Kq Tc 



(2.14) 



Equation (|2.14)l is an important content of the present 
model. This expresses how the dissipative kernel ^{t) de- 
pends on the response function ip{t) of the medium due 
to external noise e(i) [see Ea. (|2.7|l ]. Such a relation for 
the open system can be anticipated in view of the fact 
that both the dissipation and the response function cru- 
cially depend on the properties of the reservoir especially 



on its density of modes and its coupling to the system 
and the external noise source. 

To continue if we assume that e(i) is a delta correlated 
noise, i.e., (e(i)e(i')) = 2DS{t — t'), then the correlation 
function of 7r(i) is represented as 



{n{t)n{t')) - Digo^ofr-^ exp(-|t - t'\/T,) (2.15) 

where we have neglected the transient terms {t, t' > t^). 
This equation shows how the heat bath dresses the exter- 
nal noise. Though the external noise is a delta-correlated, 
the system encounters it as an Ornstcin-Uhlenbeck nose 
with the same correlation time as the internal noise 
but with an intensity depending on the couplings and 
the external noise strength. On the other hand, if the 
external noise is an Ornstein-Uhlenbeck process with 
(e(i)e(t')) = [D/t') exp(-|t - i'|/r') where D and r' are 
the strength and the correlation time of the noise respec- 
tively, the correlation function of 7r(i) is found to be 



(^(t)^(t')) 



(DgoKo) 



{r'/TaY 



— exp 

t' 



1 



1 Tc [T, 
\t-t' 



■ exp 



\t-t'\ 



(2.16) 



where we have neglected the transient terms. The dressed 
external noise 7r(i) now has a more complicated correla- 
tion function with two correlation times Tf, and r'. If 
the external noise-correlation time is much larger than 
the internal noise correlation time, i.e., r' ^ Tc, which is 
more realistic, then then the dressed noise is dominated 
by the external noise, i.e.. 



(7r(i)vr(i')) = {{Dgono)'/T'}expHt~t'\/T']. (2.17) 

On the other hand, when the external noise correla- 
tion time is smaller than the internal one, we recover 
Eq.(|233- 



III. KRAMERS EQUATION IN ENERGY SPACE 

To start with we first define the Fourier transform of 
C(i) and 7(i) as. 



7n(^) = 



dtC{t)e-Kp{-inujt), (3.1) 

dtj{t)exp(— incut). (3.2) 



In the absence of external stochastic driving force 
e{t), the fluctuation-dissipation relation {f{t)f{t')) = 
kBT"/{t — t') can be expressed in the Fourier domain as 
(now C{t - t') = (fit) fit'))) 



C';M = kBT% 

where C^(w) and jf^ are the cosine component of C„ and 
7„, respectively. Unless the explicit form of e(t) is spec- 
ified it is difficult to express the thermodynamic consis- 
tency relation H2.9|l in the Fourier domain. Without loos- 
ing generality we thus will use the general form (|3.1I3.2|) 
to derive the Fokker-Planck equation until we use explicit 
form of e(t). Conventionally low-friction regime assumes 
the relation 7 <g; w <C 1/tc, where 7 is the friction aris- 
ing due to interaction with the heat bath, evaluated in 
the Markovian limit. Tc is the correlation time of the 
noise due to heat bath and uj is the linearized system 
frequency, such relation was also considered by Kramers 
in the low-friction regime as well as for the white noise 
case. But in this paper we are concerned not only low- 
friction regime but also with non-Markovian effect due 
to bath. In this context we consider the following time 
scales in the dynamics relevant for energy diffusion in the 
non-Markovian limit l23 , 



7 < I/Tc < LU, 



(3.3) 



The separation of time scales in Eq. H3.3|) now allow us 
to write Ea. H2.4|) into the action (J) and the angle i(j)) 
co-ordinates as. 



• dx 



<j) = uj{J) 



dT-fit-T)viT)+£,it) 





dx 
dJ 



drjit - T)viT) + Cit) 



(3.4) 



(3.5) 



where v represents the velocity of the particle and for 
the deterministic part of the system's Hamiltonian, H = 

v'^/2 + Vix) we can write 



cuiJ) 



dHiJ) 
dJ ' 



(3.6) 



Since the canonical transformation (x, v) — > (J, (p) has 
been done with the deterministic part of the Hamiltonian 
it is implied that x and v can be expanded in terms of J 
and 0, 



x(J, 0) — y. Xn{J) expiin(j)) (3.7a) 

n— — 00 

oc 

v(J, 0) = ^ w„(J)exp(m0), (3.7b) 



n— — 00 



along with Xn — x!L„ and w„ = ^-n- Differentiating 
Ea. (|3.7a|l with respect to time and noting that in the 
action-angle variable space (p — toiJ) we can write 



VniJ) = inuj{J)x„iJ). 



(3.8) 



Since we are dealing with the dynamics in the one dimen- 
sion only we can choose J and (j) in such a way that we can 
make the simplification x = (1/2) X]n=-oo[^n exp(in((ii) + 
x* exp(— m(/))] for x = x* . With the choice of phase 
X = X-n (since Im(a;„) — 0), x may be further expressed 



J 



as X = J2'^=~oo ^" cos 710. Similarly using Ea. (|3.8(l we get 
^ = S^-oo "" sin ni/) for Vn = -w-n (sinccRe(u„) = 0). 
Now inserting Eas. H3.7al3.7b|) in Eas. H3.4l3.5|) we obtain 



OO CXD 



J = —i y y nXnSxp{in(j)) / dT'y{t — T)vmexp{im(j)) + i^{t) > na;„exp(i 



(3.9) 



n— — OO 771— — OO 

OO OO 






n— — OO m— — OO 



9J 



ew E 



9x^ 

a7 



exp(m0). (3.10) 



r 



In the equations H3.9|) and (|3.1U|) . the argument of the 
damping memory kernel 7 is {t — t). Now 7 decays to 
zero within the correlation time Tc- So, to deal with the 
integrals of Eas. (|3.9l3.1U|l . it is reasonable to divide the 
range of integration into two parts: (a) |i — t| ^ Tc and 
(b) i > Tc. Thus following CN |0| we can write 



m = 0[t + (t - r)] ^ 0(r) + 



dt 



{t-r), 



neglecting higher terms of Tc. It follows that 



</)(t) ~ (b{t) -{t- t)uj and v^{t) ~ v„,{t). (3.11) 

Equation H3.11|l is reasonable approximation so far as the 
of Eos. 1)3. 913. 10(1 are concerned. Within the integral, we 
therefore manipulate the behavior of (p and Vm for a time 
upto which 7(i — r) exists and also for the observational 
time at which 7 has decayed to zero. So, more specifically 
we can write for li — tI < Tc, 



00 00 



J = -i E E '^a;„Wm7m(w)exp[i(n + rn) 



n— — oc m— — OO 



+*C(^) /, fia;„ exp(i7i0). 



(3.14) 



n— — cxD 

00 oc 



dXn 

~dJ 



^{J)+ E E -^«m7m(w)exp[i( 

n— — 00 T?i— — 00 
dXn 

'dJ 



n + 771 ) 



-?(*) E ^n-expUn(? 



(3.15) 



n— — 00 



Now we are in a position to formulate the Fokker-Planck 
equation. For this we follow the method proposed by 
CN 20] which is based on Kramers-Moyal expansion of 
the transition probability that connects the probability 
distribution function P(J,(j),t) at time t with that of 
P{J, (pjt + t) at a later time t + t for small r, given 
that we know the moments of the distribution. The time 
evolution of the probability distribution P{J,(j),t) is de- 
termined by the equation. 



dTj(t — T)Vm{T) exp[i7Tt0(T)] 



Vm{t) exp[im<f>{t)] / dTj{t — t) exp[—im{t — t)llj] 



(3.12) 



and for t ^ Tc, using Ea. ((3.2|l we have 



dTj{t — T)vmiT) exp[i?Ti0(r)] 



Vm{t) exp[i7TT.(/)(t)]7„(tj). 



(3.13) 



Putting Ea. H3.13() which takes into account the observa- 
tional time scale, in Eas. (|3.9|) and H3.10|l we get 



dP 

'dt 



lim 

T-+0+ 



_. OC 

7E 



(-1)" 



E 



d 



n=l 

k 



m,k—0-^{m+k—n) 



d_ 
dJ 



X — mAJtr{A^,r))p} 



(3.16) 



where AJt = AJt(r) = J{t + t) - J{t) and A0t = 
A(/>t(T) = (/)(t-|-r) — (/)(t). At this juncture it is worth re- 
calling that r is the coarse-grained time scale over which 
the probability distribution function evolves, whereas t^ 
is the correlation time, which due to low damping is much 
smaller than t. The low value of 7 prompts us to take 
1/7 as the largest time scale for the entire problem. How- 
ever, the reciprocal of the frequency of oscillation, i.e., 
l/w, is the smallest time scale. Our task is now to cal- 
culate the moments of the form (((AJ()™(A0()'')) where 



((• • ■ )) means that we have taken the two averages inde- 
pendently. 

To evaluate the moments we make use of the following 
standard procedure [23, uM 



AJtir) 



dsJ[Jit + s),(t)(t + s),t + s] (3.17) 



A</)t(r) = / ds(j)[J{t + s),(l>{t + s),t + s] (3.18) 
Jo 

where J and (p are given by Eas. (|3.14(l and H3.15|l . re- 
spectively. 

The non-Mar kovian nature, i.e., r^ is finite but t^ < t, 
of the present problem allow us to consider all orders 
of r in Ea. H3.16|) . But, since dP/dt is evaluated in the 
limit T — > 0+, terms linear in t are taken while all the 
higher powers are neglected. We then recast Eas. H3.14|) 
and H3.15|) in the following form 

oo oo 

i = - 2J z2i -S„m(J)exp[i(n + r7i)(/)] 

n— — oo m— — oo 
oo 

H{t) Y. ^n(J)exp(z7i0) (3.19) 

n— — oo 

OO OO 

(j) = w(J)+ ^ ^ C„™(J)exp[i(n-|-m)^] 

n— — OO 7n— — oo 

oo 

-?W Y. Mn(J)exp(m</.) (3.20) 

where 



n— — oo 



cr„(J) = inXn{J), 
dXniJ) 



^J■n{J) 



dJ 



(3.21) 
(3.22) 



Bnm{J) = inXn{J)v„i{J)^^[uj{J)] 
dXn{J) 



^nni \'^ ) 



dJ 



(3.23) 

Vm{J)lmHJ)]. (3.24) 



Finally the moments can be calculated using the stan- 
dard iterative process prescribed by CN |23 and they are 
of the following form 



((AJt(r)» 



-2r£n2[c.|x„p7^(c.) 



-— {|a:„pa;i(u;)} 



(3.25) 



((A0t(r))> 



UJT + T y n 

n=l 



■E^ 



^%{^) 



d\Xr^ 



dJ 



d (P,s(.d]xnl_ 

'dj[^-^''>^r 



(3.26) 



{{[AMr)f))^ArY,n'\x„\'d^^{co), (3.27) 



(([A0,(r)]2))=4r£r 
((AJ,(r)A0,(r)))=O, 



UXfi 



dJ 



C;j(c.), (3.28) 
(3.29) 



where 



In = 



Cc 
r, 

/o 



ai 



dt'-f{t) cos{nujt), 


(3.30a) 


dtj{t) sin(na;t). 


(3.30b) 


dtC{t)cos{nujt), 


(3.30c) 


dtC{t)sm{nujt). 


(3.30d) 



Also 



%iLo) = %{Lo)-tr„W), (3.30e) 

dnicj) = C^^{uj) - iC^M- (3.30f) 

In the absence of the external noise e(t), C'„(aj) reduces to 
CJlo) = kBT%{Lu) for which Eas. l|X^ - (|T^ becomes 



(AJ,(r)> = -2r^n2 L _ fc^T— ) (K^T^i), 



(A(/>t(T)) = wr -}- r N^ n luj — ksT 



dJ 



d\x„ 



dJ 



-In 



([AJt(r)n=4rfcBr£ 

oo 

([A0t(r)]2)=4TfcBT^: 

n=l 

{AJt{T)AMT))=0- 



n^\Xn?ln. 



dXr. 



dJ 



Note that in the above unnumbered equations there is 
only one averaging, (• • • ) instead of two averaging, ((•••)) 
used in this article. This is due to the fact that in the 
present model we make an extra averaging over the ex- 
ternal noise processes in addition to the usual thermal 
averaging procedure. 

Inserting Eas. (|3.25|l - H3.29|l in Ea. (|3.16|l and neglecting 
terms with n > 2 we obtain the Fokker-Planck equation 
for P{J, (/>, t) as 



dP{JA,t) 
dt 



d_ 
dJ 



^W S?MiL-w p 






dP 



(3.31) 
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where 
















e(J) = 


oo 




(3.32) 


r(j). 


oo 

n— 1 


dXn 
dJ 


2 




(3.33) 


fi(J) = 


oo 
n=l 


n 




dJ 


[" dJ )\ 





(3.34) 

For a distribution function that is initially (t = 0) inde- 
pendent of (j) the diffusion equation in action space be- 
comes 



frequency of the dynamical system, i^iE). Although the 
expression for diffusion coefficient (|3.4(J|) looks bit com- 
plicated and formal due to the appearance of the Fourier 
coefficients x„ in the summation, it is possible to read 
the various terms in D{E) in the following way. D{E) 
is essentially an approximate product of three terms, 
A(wo), j^ dt-f{t)cos[nv{E)t], and J^(^) I],?=i "^l^^"^' 
where the n dependence of the latter two contributions 
have been separated out for interpretation. The first 
term is an analogue of ksT for the open system, the 
integral is the Fourier transform of the memory kernel, 
while the sum can be shown to be equal to J (see Ap- 
pendix D of CN [23|), the action variable. For a system 
only coupled to a heat bath, i.e., for no external driving, 
D{E) reduces to the expression derived by CN J23|. 



dP{J,t) _ d 



dt 



where 



e{J)\K— + u:{J)\P 



A = A(wo) 



%(wo) ■ 



(3.35) 



(3.36) 



Here wq is the linearized frequency and A plays the typi- 
cal role of ksT which for e{t) = becomes equal to ksT . 
Now by virtue of Eq.(|S3I) w( J) = dH/dJ = dE/dJ. Ex- 
pressing 



IV. ENERGY DIFFUSION CONTROLLED 
RATE OF ESCAPE 

The classical treatment of memory effects in the energy 
diffusion controlled escape is now well documented in 
the literature [23, |33, IMI . To address the corresponding 
problem for the open system we first rewrite the Kramers 
equation l|3.39|l in the form of a continuity equation 



dP{E,t) dJE 

dt BE 



(4.1) 



uiJ) = iy{E), 



where je is the stationary flux along the energy coordi- 
(3.37) nate and is given by 



we have 



^-<-)^- 



(3.38) 



With this transformation, for an external noise driven 
bath, the Kramers equation for energy diffusion 
[Ea. H3.35|) ] becomes 



dP{E,t) _ d 
dt ~ 'dE 



D{E)[^ + \].{E)P{E,t) 



(3.39) 



with the following diffusion coefficient 



D{E) ^ u{E)2K{ujQ)y^n^\xr,\^ / 



dt-f{t)coii[7w{E)t\. 



(3.40) 

Eq. ^3.yy\) is the first key result of the present article. The 
equation is valid for arbitrary temperature and noise cor- 
relation. The prime quantities that determine the equa- 
tion for energy diffusion H3.39|l are the diffusion coeffi- 
cient D; the open system analogue of ksT, A; and the 



JE = -D{E) 



d_ l_ 
dE^A 



HE)Pst{E), (4.2) 



with Pst being the stationary probability distribution. 
For zero current condition, we have the stationary distri- 
bution, Pst at the source well 



Pst{E) = ^j^^cxp{-E/A) 



(4.3) 



where N is the normalization constant. Here it is impor- 
tant to mention that for e(i) = 0, one has Pst = Peg- We 
now define the rate of escape k as fiux over population 

M 



k = JE/ria (4.4) 

where ria is the total population at the source well, 



n„ = / P{E)dE. 



(4.5) 



Here Ei, is the value of the activation barrier. Follow- 
ing Biittiker, Harris, and Landauer (BHL) [3y| we use a 
Kramers like ansatz 



Setting r]{E) — r]{Eb) exp[s(£' — Eb)/A] and putting this 
into Eq. (|4.7I) and comparing this with the right hand side 
of Eq. (|4.8() we have 



P{E)^r^{E)pst{E) 



to arrive at 



(4.6) viEb) = 1/(1 - s) for s < 0. 

Thus escape rate k can be written as 



(4.13) 



JE ^ -D{EME)pst{E) 



dyjE) 
dE ' 



(4.7) 



Integrating the above expression from E — Ei c:^ K io 
E = Eb, one derives an expression for energy independent 
current je (with E ^ E^) as 



k ^ JE 



f]{E)pst{E)dE 



(4.14) 



Making use of Ea. H4.13|) in Ea. l|4.8l) and the resulting 
expression for je in Ea. (|4.14|l we obtain 



jE 



77(A) - r^iEb) 



L 



Eb 



dE 



A D{E)u(E)pst(E) 



= [\-7^{Ei,)]D{Et,)^^eM-Eb/K), (4.8) 

where we have used the boundary condition 77(A) ~ 1. 

Following the original reasoning by BHL we now allow 
an outflow jout from each energy range E to E + dE, 
with each E satisfying the condition E ^ Eb- Then we 
can write 



djout = ai^iE)rj{E)p,tiE)dE, (4.9) 

which is compensated by a divergence in the vertical flow 



dJE 
dE 



ay{E)f^{E)pst{E). 



(4.10) 



Here a is a parameter that has been set approximately 
equal to one by BHL, though in g eneral the parameter a 
is not always equal to one p, l37l | . Inserting the expres- 
sion for nonequilibriuni current fEa. l|4.7|l ) in the above 
expression we obtain an ordinary differential equation for 
V{E) 



d'^rj 



D{E)^2 



dD D{E) 

'dE A~ 



^-ari{E) = Q. (4.11) 



fc = 



1-s 



j^' ri{E)p,t{E)dE 

{N-^/A)D{Eb)cxp{-Eb/A) 



(4.15) 



For the dynamics at the bottom we have 77 — > 1. For luq 
being the frequency at the bottom of the source well we 
now calculate the total population of the source well. 



00 /-OO 



Tia = / / Pst{E)dxdp 

-OQ J —00 

= iV-i(27rA/wo)- 



(4.16) 



So, for the external noise driven heat bath the non- 
Markovian rate of escape from a metastable well in the 
low friction regime is given by 



{1 + {AaK^/D{EbW^ - 1 
{1 + (4aA2/^(^fc))}i/2 + 1 

xwoexp(-£;{,/A). 



D{Eb) 

A2 



(4.17) 



Eq. i4.17[) is the second key result of present paper. 



V. SPECIFIC EXAMPLE: HEAT BATH DRIVEN 
BY EXTERNAL COLOR NOISE 

As a specific example, we consider that the heat bath 
is modulated externally by a colored noise e(i) with noise 
correlation 



Within small energy range above Eb one can as- 
sume essentially a constant diffusion coefficient, i.e., 
dD{E)/dE\E~Ei = for ^ ^ Eb. Now substituting a 
trial solution of the form ri{E) ~ '^exp(si?/A) for s < 0, 
in Ea. H4.11|) we have 



4aA2 
D{Eb) 



1/2 



(4.12) 



(e(i)e(i')) = — exp 



\t-t'\ 



(5.1) 



where D^ and Te are the strength and the correlation 
time of the external noise, respectively. In addition to 
that we also consider the internal noise /(i) to be white. 
The effective Gaussian Ornstein-Uhlenbeck noise ^(i) = 
f{t) + 7r(i) will have an intensity Dji and a correlation 
time Tfl given by [l5j 



Df 



{mm)dt, 



^R = ^ /°°i(e(i)^(0))dt. 

^B Jo 



(5.2) 



(5.3) 



Following the above definitions and using Ea. H2.17|l we 
have 



Dr = gliksT + Denl) and tr 



Dr 



(5.4) 



It is important to mention here that since we are treat- 
ing the internal noise processes to be a delta correlated 
one(rc —^ 0), Tc does not appear explicitly in the expres- 
sion of Dfi and tr. With this the eflFective noise ^{t) 
becomes a colored noise and its correlation is given by 



((e(Oe(t')» = — exp 

TR 



\t~t'\ 
TR 



(5.5) 



To study the dynamics we consider a model cubic po- 
tential of the form V{x) — Ax^ — Bx^ where A and B are 
two constant parameters with ^ > and B > Q. The dif- 
fusion coefficient D(Ei,) in the internal white noise limit 
reduces to 



DiEb) = 5oA(c^o) J 
where the action, J is represented as [2C 



(5.6) 



J^2iy{Eb)Y,ri^\^n\ 



(5.7) 



and can be calculated using the following standard form 



TT 



V dx 



(5.8) 



where Xi and X2 are the two turning points of oscillation 
for which v is equal to zero and they both corresponds 
to total system energy E. In principle, they are the first 
two roots (in ascending order of magnitude) of the cubic 
equation 



-Ax^ + Bx^ + E = 0. 



(5.9) 



For an external color noise driven heat bath A(ujo) (see 
Ea. (|3.36|l ) reduces to 



We then numerically solve the Langevin equation l|2.4|) 
using the second order stochastic Heun algorithm ,38. .39.] . 
To ensure the stability of our simulation we have used 
a small time step At = 0.001 with At/rji < 1. The 
numerical rate has been defined as the inverse of the mean 
first passage time [40, |4l| . The mean first passage time 
have been calculated by averaging over 5000 trajectories. 
The value of the other parameters used are given in the 
caption of Figs.(l) and (2). 

In his dynamical theory of chemical reactions Kramers 
identified two distinct regimes of stationary nonequilib- 
rium states in terms of dissipation constant (7). The es- 
sential result of Kramers theory is that the rate varies lin- 
early in weak dissipation regime (characterized of diffu- 
sion of energy) and inversely in the intermediate to strong 
damping regime (spatial diffusion limited regime). That 
is, in between the two regimes the rate constant as a func- 
tion of dissipation constant exhibits a bell-shaped curve 
known as Kramers turnover _4, _5j . In the traditional sys- 
tem reservoir model the dissipation and the fluctuation, 
both originating from a common source, the reservoir, are 
connected through the fluctuation-dissipation relation. A 
typical signature of this relation can be seen through the 
turnover phenomenon in Kramers' dynamics. Whereas 
for a thermodynamic open system where the heat bath 
is modulated by an external noise, both the dissipation 
and the response function depend on the properties of 
the reservoir, mainly on its density of modes and its 
coupling to the system and the external noise source. 
By virtue of this connection between the dissipation and 
the external noise source, Ea. (|2.9|l plays the typical role 
of the thermodynamic consistency relation, an analogue 
of the fluctuation-dissipation relation for thermodynamic 
closed system, for which one can expect turnover like fea- 
ture in Kramers' dynamics (for the open system). So, for 
the external color noise driven bath we first wanted to 
check whether Kramers' turnover feature can be restored 
from our model. In Fig. 1 we have plotted the rate con- 
stant, k obtained from Langevin simulation, for a wide 
range of damping constant, Qq for different values of ex- 
ternal noise correlation time, Tg. The figure shows usual 
Kramers' turnover of the rate constant with variation of 
the damping constant. The shift of the maxima occurs as 
the external noise correlation time varies, a typical effect 
of the bath modulation. 

Next we compared the theoretical result (|4.17|) with the 
numerical simulation data. In Fig. 2 we have plotted the 
rate constant, k against the damping constant, g^ in the 
weak damping domain (0.001 ^ g^ ^ 0.01) for different 
values of the external noise correlation time, Tg. What 
we observe is that the agreement between the theoretical 
prediction and numerical simulation is quite satisfactory. 



VI. CONCLUSION 



A(wo) = 



ksT + D,4 



1 -I- w^r- 



2.^2 



(5.10) 



Based on a simple system-reservoir Hamiltonian ap- 
proach, we have studied the behavior of a subsystem cou- 
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FIG. 1: Turnover phenomenon for external color noise driven 
bath. Parameters used are ksT = 0.1, De = 1.0, Kq = 5.0, 
a = 1.0, A — 0.5 and Et = 5.0 (scale arbitrary). 
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FIG. 2: Barrier crossing rate in the low friction regime 
(0.001 ^ go ^ 0.01), a comparison between theoretical pre- 
diction, Eg 14. 171 (solid lines) and Langevin simulation. Pa- 
rameters used are ksT — 0.1, De = 1-0, kq — 5.0, a — 1.0, 
A = 0.5 and Et — 5.0 (scale arbitrary). 



state Kramers' escape rate from a metastable well in the 
low friction regime. The main conclusions of the present 
work are the following: 

(i) Since the reservoir is driven by the external noise 
and the dissipative properties of the system depend on 
the reservoir, we have established a simple relation be- 
tween the dissipation and the response function of the 
medium due to external noise. This relation is impor- 
tant for identifying the effective temperature of the heat 
bath characterizing the stationary state of the thernio- 
dynamically open system. 

(ii) We then followed the dynamics of the open sys- 
tem in the energy space and derived the corresponding 
Fokker-Planck equation with diffusion coefficient contain- 
ing the effective temperature like quantity which is an 
open system analogue of ksT. Following the standard 
approach we then derived the generalized non-Markovian 
Kramers' escape rate from a metastable well in the en- 
ergy diffusion regime. 

(iii) From the point of view of the realistic situation 
we considered the special case where the internal noise 
is white and the external noise is colored and have cal- 
culated the escape rate for a model cubic potential. We 
have shown that the theoretical prediction agrees reason- 
ably well with numerical simulation. In addition to that 
we have also shown that our model recovers the turnover 
feature of the Kramers' dynamics when the external noise 
modulates the reservoir. 

(iv) Finally, as shown in the present work one can eas- 
ily tune the external noise parameters from outside which 
can be used to study the effect of several kinds of noise 
properties, e.g., long tail gaussian noise j3], in Kramers' 
dynamics. Another suitable candidate for studying the 
escape rate dynamics can be irreversibly driven environ- 
ments |42L 14!^ . In our future communication we would 
like to pursue such theoretical analysis. 
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